Tunable electronic transport and unidirectional quantum wires in graphene subjected 

to electric and magnetic fields 
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Magnetic barriers in graphene are not easily tunable. Here we show that the application of both 
electric and magnetic fields provides tunable and far more controllable electronic states in graphene. 
In particular, a one-dimensional channel (quantum wire) can be created, which supports localized 
electron-hole states with parameters tunable by the electric field. Such quantum wire offers peculiar 
conducting properties, like unidirectional conductivity and robustness to disorder. Two separate 
quantum wires comprise a waveguide with two types of eigenmodes: one type is similar to traditional 
waveguides, the other type is formed by coupled surface waves propagating along the boundaries of 
the waveguide. 
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PACS numbers: 71.10.Fd, 73.20.Fz, 73.23.-b 



I. INTRODUCTION 

The realization of stable single-layer carbon crystals 
(graphene) triggered an explosion of interest in this ma- 
terial because of its unique electronic properties (for re- 
views see, e.g., [J-Q), making it a promising candidate 
for designing one-chip nanoelectronic devices (e.g., 0,0- 
7]). However, Klein tunneling Q hinders the applica- 
tion of traditional methods of current control (e.g., on- 
off switching, changing the current direction, etc.) by 
tuning the voltage between various elements of a device 
. This effect also complicates the creation of localized 
electron-hole states in graphene. 

As it was shown recently [l(J E| , this difficulty in cre- 
ating localized states can be overcome by using an inho- 
mogeneous magnetic field. An alternative way (which is 
mathematically analogous to the previous one), involv- 
ing graphene sheet deformations, was proposed in [T^ - 
|X4| ] . The regions with either inhomogeneous magnetic 
field or strains can act as non-trans par ent barriers and 
angle-resolved charge carriers filters [l5| and form long- 
lived [l6[ or stationary localized states El, G3, 13 ■ 
These 1-D states appear as "surface" waves p ropaga ting 
along the inhomo gen eity of the magnetic field [19l - l24j , or 
along the strains [12| . 

The charge-confining and guiding capabilities of mag- 
netic barriers [l(| [H| and graphene strains [HI open 
up certain possibilities for potential future applications. 
However when it comes to designing fast-tunable elec- 
tronic devices (switches, filters, etc.) a difficulty emerges. 
The problem is that most of the existing magnetic barrier 
technologies usually imply the deposition (on top or be- 
neath the graphene sheet) of a fixed pattern of magnetic 
material which reproduces the desired magnetic field dis- 
tribution in the sample. Any change of parameters means 
in fact building a new setup and creates formidable (if 
surmountable) obstacles for harnessing magnetic barri- 



ers as elements of fast-acting electronic devices. In other 
words, magnetic barriers are not easily tunable. 

Here we suggest an efficient way around this prob- 
lem by simultaneously employing both inhomogeneous 
magnetic and electric fields: i.e., combined electric and 
magnetic barriers. The proper combination of these two 
allows a better control of the transport properties of 
graphene by tuning the electric potential, with the pa- 
rameters of the magnetic field remaining intact. Depend- 
ing on the voltage, this barrier can be either semitrans- 
parent or opaque. 

The combined electrostatic and magnetic barrier pos- 
sesses a unique feature that makes it different from other 
types of barriers. Graphene with mutually perpendic- 
ular electric and magnetic fields supports states which 
are localized near the barrier. These current-carrying 
states (surface waves) correspond to quasiparticles mov- 
ing along the barrier just in one direction. The direc- 
tion and the value of the quasiparticle velocity is easily 
controlled by the electrostatic potential. These states 
correspond to the classical drift of charged particles in 
crossed electric and magnetic fields and exist if and 
only if the drift velocity is smaller than the Fermi ve- 
locity. The absence of counter-propagating states pre- 
vents the backscattering induced by either irregularities 
in graphene [25l427| or by the fluctuations of the mag- 
netic field. For more about the effects of disorder on the 
electronic properties of graphene see Chapter IV in [1] 
and references therein. 

While one barrier forms a wire, two such barriers make 
up a waveguide. This waveguide has a set of "ordinary" 
waves and another set of "extraordinary" waves. The 
ordinary waves are characterized by the quantization of 
their transverse wave numbers, while the extraordinary 
waves are formed by two coupled surface waves propagat- 
ing along the waveguide walls (barriers). Depending on 
the barrier parameters, the extraordinary modes can be 
either bidirectional or unidirectional. There is an energy 
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gap where only extraordinary modes exist. Decreasing 
the spacing between the barriers broadens this gap. The 
extraordinary modes are also stable against backscatter- 
ing. 

The unique, easily controllable, and tunable features 
of the combined barriers and waveguides are promising 
toward creating new graphene-based electronic devices. 

This paper is organized as follows. In Sec. II, the mod- 
els for the electromagnetic barriers and auxiliary con- 
structions are described. In Sec. Ill, the properties of a 
single electromagnetic barrier are studied. The spectrum 
and eigenfunction of the waveguide formed by two barri- 
ers are presented in Sec. IV. Concluding remarks are in 
Sec. V. 



II. BASIC FORMALISM 
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A. Step-like barrier. 

The low-energy excitations in single-layer graphene in 
the presence of a perpendicular to the layer magnetic field 
H z (x) = dAy/dx and an in-plane electric field E x {x) = 
—dV(x)/dx (which are constant along the y-direction) 
are described by the Dirac equations: 



i - + - + iC A(x^ 
dx dy ch ' 
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dy 



ch 



A(x) 



tpA 



[£ - eV{x)]il> A {l) 



where A(x) is the y-component of the vector poten- 
tial, V(x) is the scalar potential, £ is the energy of 
the quasiparticle, and ip is the two-component spinor 

Let us consider two homogeneous graphene domains 
subjected to different constant scalar (V% and V2) and 
vector (A y i = A\ and A y2 = A 2 ) potentials, and assume 
that the domains are connected by an inhomogeneous 
transition region where the potentials vary smoothly be- 
tween the constant values (see Fig. Q]). When the width 
I of this region is large compared to the graphene lattice 
spacing a, and small compared to the Fermi wavelength 
Xf (a <C i <C Ajr), the potentials can be replaced by step- 
like functions, as shown by the red dashed line in Fig. [T] 
This widely used approximation simplifies the problem 
considerably. 

Due to translational invariance along the y-direction, 
the solutions of Eqs. (JXJ) can be presented in the form 
il>A,B(x-,y) = exp(ik y y)ty a.b(%)- In dimensionless vari- 
ables £ = x/L, 



K± 



kyL^ s 



£L/hvF, 



u = eVL/hvF, A—(eL/ch)A, 
where L is a characteristic spatial scale (magnetic length 



FIG. 1: (color online) Schematic diagram of the inhomo- 
geneity region between two otherwise-homogeneous graphene 
sheets. The function P(x) can be either the electric scalar 
V(x) or magnetic A(x) potential (hence the choice of P). 
Here, P(x -> -00) ->• P u and P(x -> 00) -> P 2 . The 
dashed line is the step-like barrier approximation. Note that 
A = A y y = A{x)y. 



Ih = ych/eH, for instance), Eqs. ([1} have the form: 



— - kx - A ] 4>a = *(e - u)4>b, 



— + k_l + A) tp B = i(e - u)ip A . 



(2) 



B. Homogeneous graphene. 

In homogeneous graphene (A — constant, u = 
constant), $ a,b(£) ~ ex P(* K ||£): an d the wave vector 
components ki\ and k± are related by the dispersion re- 
lation: 



k\ + (k ± + Af = (e- uf 



(3) 



This equation is valid for both propagating (Im/tii = 0) 
and non-propagating (evanescent) waves (Re/C|| =0). In 
the («|| , Kj_)-plane, the wave vectors k of the propagating 
waves lie on a circle centered at the point kj_ = —A, with 
radius p = \e—u\ (Fig. [5]). For given k± and e, there are 
two solutions of the dispersion equation ([3]) with positive 
(blue solid arrow in Fig. [5]) and negative (blue dashed 
arrow in Fig. [2} values of 



C. Probability current. 

The solutions of Eq. @ have the form: 

^AM^^B^ + ^A^' 1 ^- 



(4) 
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FIG. 2: (color online) The wave vectors k (blue thick arrows) 
of propagating waves in homogeneous graphene lie on a circle, 
described by the dispersion equation ([3|. The normalized 
coefficients pC^ + ' (red solid arrow) and pC' - ' (red dashed 
arrow) are shown when sgn(e — u) = +1. The current with 
positive J x component is directed either along the solid red 
arrow when sgn(e — u) = +1, or along the green dotted arrow 
when sgn(e — u) = —I. 



The amplitudes are connected by the relation: 



where 



and 



C(±) = 



i(K± + A) ± Kll 



c<-> = -1/CW. 



(5) 
(6) 
(7) 



Eq. ([5]) allows the consideration of the spinor component 
A only (thereafter \E'). When the wave vector k lies on 
the circle (propagating waves), the following representa- 
tion of the coefficients ' is valid: 



C* (+) = C = sgn(e - u)e^, C<-> = -C* 



where 



Sill (y9 : 
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COS If = 



(8) 



(9) 



The angle </? is connected with the direction of the proba- 
bility current density, J [1(3, EH ■ Indeed, the probability 
current density J — {i]j\cr\ip) can be presented in the 
form: 

J x + iJ v = 2ip* A ip B , (10) 

Using Eq. (|T0|) . the current densities that corre- 

spond to pure (+) and (— ) states, can be written as: 

^(+) 



J (+) = 2sgn{e-u)e ltp 
= — 2 sgn(e — u)e 



(-) 



It follows from Eq. (ITT1) that the current with positive 
component J x > is described by ?/^ + ) when sgn(e — u) = 
+1 (red arrow in Fig. [TJ, and by when sgn(e — u) = 
— 1 (green arrow in Fig. [TJ. 



III. COMBINED ELECTRIC AND MAGNETIC 
BARRIER 

A. Semi-transparent barrier. 

The continuity condition for the spinor components 
^ a,b on the interface between two homogeneous domains 
(domain 1 and domain 2) can be written in the form: 



,v4" hT 



(+) ./.(-)\T 



where M is the transfer matrix: 
1 



M = 



C 2 + C^ 1 



C 2 1 + Ci C 2 

C*2 — Cl C*2 



cr 1 
cr 1 



(12) 



(13) 



with 



The coefficients are defined by Eqs. ©-([8]), 
potentials u\^% and A\.2, correspondingly. 

To describe the transport properties at the barrier, we 
now introduce a graphic representation, shown in Fig. |3l 
which provides a better understanding of the dispersion 
relations. In Fig. [3l points on the circles represent solu- 
tions of the dispersion equations in domains 1 and 2, and 
correspond to propagating waves (waves with real ftii) in 
these domains. 

Let us now consider the "refraction law" which relates 
the directions of the incident and refracted waves [states 
with the same sign of the current component J x in both 
domains (for definiteness, J x > 0)]. This "refraction law" 
can be easily derived with the help of Fig. [3] as follows. 
Due to translation invariance along the y-axis, the wave 
vector component k± has the same value in both me- 
dia. The wave vectors At of the propagating waves take 
on values n\ and Ki in the first and second domains. 
The wave vectors K\ and lie at the intersection of 
the corresponding circles and the line k±_ = constant 
(thin horizontal dashed black line in Fig. [3]). There 
are two such intersections for each of the two circles. 
The physically meaningfull intersections (i.e., the inter- 
sections that present the solutions with positive current 
density component J x > 0) are defined by the sign of 
(e — Mi, 2)- The currents J are directed cither along the 
radius vector (i.e., the vectors from the centers of both 
circles to the intersections, in Fig. of the intersection 
or in the opposite direction, depending on the sign of: 
(e — Mi. 2). For example, when sgn(e — u±) = +1, the inci- 
dent current J is directed along the red arrow in Fig. [3] 
The refracted current is directed along the blue arrow 
when sgn(e — u 2 ) = +1, or along the green arrow when 
sgn(e - u 2 ) = -1. 

The following relation defines the connection between 
the incident, 6\, and refracted, 62, angles: 



(11) 



K± 



-Ai - (e — ui) sin 0i = — A2 - {s — U%) sin# 2 . (14) 
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FIG. 3: (color online) Here k± and K\\ are the transverse and 
longitudinal components of the wave vector n. The circles 
represent solutions of the dispersion relations e(fc) in the do- 
mains 1 and 2, and correspond to propagating waves. The 
top circle has a radius p2 = |e — ita], while the bottom circle 
has a radius p2 = \s — Ua|. The incident current is directed 
along the bottom-right red arrow. The refracted current is 
directed along the blue arrow (the only arrow pointing South- 
East, towards 5 o'clock), when sgn(e — U2) = +1, and along 
the green solid arrow (pointing North-East, 1 o'clock) when 
sgn(e — U2) = — 1. The latter one is opposite to the dashed 
(South- West, 7 o'clock) green one, which intersects the black 
horizontal dashed line. 



When, forinstance, |e — u\\ = \e — ua| = | £ — u |) the 
refraction law reads: 



sgn(e — u\) sin#i + sgn(e — u%) sin #2 



(15) 



The refraction law presented in [13, HH is a particular 
case {u\ = u 2 — 0) of this relation. 

It is possible to show that the transmission coefficient 
T = J x \/J x 2 is denned by the following expression: 



T = T(w 2 ,wi) 
where 



2 V / r^ 2 yr 



1 + y/l — w\ y/\ — w\ - W 2 W\ 



(16) 



w 2 ,i = 



A 2 ,i 



e - u 2 ,i 
It follows from Eq. JTB} that 



T\ 



1. 



(17) 



i.e., there is an angle of incidence, non- normal in the gen- 
eral case, for which the interface is totally transparent. 
In other words, the barrier changes the direction of inci- 
dence, at which Klein tunneling occurs [28j . 



B. Non-transparent barrier: bound states. 

When the difference between the vector potentials is 
large enough, 



\Ai - A 2 \ >\e- m\ 



u 2 \ 



(18) 



the circles in Fig. [3] do not cross, and no propagating 
wave can penetrate through the interface. Therefore, the 
transmission coefficient is equal to zero, T = 0, and the 
barrier is nontransparent (reflecting wall). However, it 
can support a wave, which propagates along the line sep- 
arating two domains with the amplitude exponentially 
decaying in the transverse directions. This mode is a 
ID analog of two-dimensional surface waves. In order to 
simplify terminology, in what follows, we call it a surface 
wave. 

Let us now determine the properties and existence con- 
ditions of these surface waves. When the wave vector 
components K\\j (J — 1,2) in both media are imaginary, 
«ll j = i |«|| j| the wave functions ^1,2 have the form: 



HL+hiil£ 



(19) 



where k\\ j — \J{e — Uj) 2 — /tj_ = ~ ( £ ~ u j) 2 ■ The 



wave is now localized near the interface 



',(-) 



the condition 
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0, i.e. 



Mil 



4 +) 



(+) 



dC 2 + 1 
Cf + 1 







£ = when 
0. Therefore, 



(20) 



is the dispersion relation of the surface waves. A solution 
e(«x) of this equation exists if and only if 



til — u 2 



Ai - A 2 



< 1 



and has the form: 



s(ka 



uiA 2 - u 2 Ai 
A 2 - Ai 



Ul - u 2 

A2-A1 



(21) 



(22) 



Equation (|22)) describes a surface wave propagating along 
the line £ = 0, with the group velocity v g ~ de/dn± = 
(ui -u 2 )/(A 2 - Ai). 

Note that the inequality (|2"Tj) in the dimensional vari- 
ables takes the form: 



c 

v F 



AV 



AA 



< 1, 



(23) 



where AV and AA are, respectively, the differences 
between the scalar and vector potentials in the two 
graphene domains. From Fig. [T] one can see that the elec- 
tric field in the inhomogeneous region is E x ~ —AV/i, 
and the magnetic field is H z ~ AA/£. Therefore, the 
inequality (|23j) can be written as: 



Vd = c 



Ex 



< Vf- 



(24) 
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Thus, the dimensional group velocity v g — h~ 1 d£/dk y 
of the surface wave, coincides with the drift velocity 
Vd = cE /H of a charged particle in crossed electric and 
magnetic fields. 

If, for example, the characteristic width, £, of the 
barrier (see Fig. [TJ is £ = lOnm, and the magnetic 
field H z = 1 T, the condition for the surface wave to 
exist, Vd < vf [Eq. (|24[) ]. is satisfied when the potential 
difference across the barrier (between two graphene do- 
mains in Fig. [T]) is AV = lOmV or less. These numbers 
are quite feasible. Magnetic barriers with amplitudes of 
up 1 T have been created experimentally by depositing 
ferromagnetic films on top of a graphene sheet |29l |30| . 
Patterened stripes down to lOnm can be realized using 
nanolitography (see (3ll | for review). 

Note that the bound state near the <5-function magnetic 
barrier [H z ~ 5(x)] described in is the particular case 
Ui = U2 (E x = 0) of this surface wave with zero group 
velocity. 



C. Regions in the (ftx,e)-plane where surface waves 

exist 

The regions in the (k±, e)-plane where the surface 
waves exist can be defined using the graphic construc- 
tion shown in Fig. 2) The Dirac point (— Ai,u\) in the 
graphene domain 1 generates a division of the plane into 
four sectors. Two of them (yellow sectors in Fig. @Jt) 
correspond to propagating waves with real rey, and the 
other two sectors correspond to non-propagating waves 
with imaginary km (white regions in Fig. H£i). A simi- 
lar division of the plane is generated by the Dirac point 
{—A2,U2) in the domain 2. The blue sectors in Fig. 0k 
correspond to propagating waves in this domain. 

Green regions in Fig. 0^ show the overlap of the 
propagating- wave sectors, i.e., waves whose parameters 
k_l and e lie in these overlapping green regions can simul- 
taneously propagate in both domains. For these waves 
the interface between the graphene domains acts as a 
semitransparent barrier. 

Yellow (blue) sector corresponds to the waves that can 
propagate only in domain 1 (only in domain 2), and the 
interface between the domains acts as a non-transparent 
(reflecting) barrier. Parameters «j_ and e of the non- 
propagating surface waves (i.e., waves that are evanes- 
centin the both domains) belong to the white regions 
in Fig. 5^,. The red solid line that connects the two 
Dirac points is the dispersion curve e(n±) described by 
Eq. (j22j) . It follows from Eq. (l22j) that the surface waves, 
if exist, propagate in only one direction. In other words, 
such barrier constitutes a unidirectional quantum wire. 

When the inequality (f2"Tj) is not satisfied, the divi- 
sion of the (k_l, e)-plane has different structure, shown in 
Fig. HJd. Here there is no white region between the Dirac 
points, i. e., the surface waves are absentg in this case. 
Another important difference between these two cases is 
the following. When the inequality (f2"Tj) holds, there is 



(a) 











(b) £ 











FIG. 4: (color online) The yellow and blue sectors corre- 
spond to propagating waves in domains 1 and 2, respectively. 
The overlapping (green) regions correspond to waves that can 
propagate both in domain 1 and domain 2. (a) the inequality 
Eq. (|21|) is satisfied. The barrier is opaque for all the waves 
whose energy lies in the energy gap between the green re- 
gions. The dispersion e(ki) of the surface wave is shown by 
the straight red line joining both Dirac points, (b) the in- 
equality Eq. (|21|) is not satisfied and there is no surface wave. 
The barrier is semitransparent, i.e., regardless of the energy, 
there are always the angles of incidence when the waves pen- 
etrate through the barrier. 

an energy gap, in which there are no waves penetrat- 
ing through the barrier, and the barrier is opaque for all 
angles of incidence (all k±). Figure |4ji allows easy deter- 
mination of this range of energies. The absence of waves 
penetrating through the barrier means that the horizon- 
tal line e = const does not cross any green regions. It 
is readily seen that this condition is satisfied when the 
energy lies in the following range: 

i(wi + u 2 - \Ai -A 2 \) <e< + u 2 + \Ax -A 2 \). 

When the inequality (f2"Tj) is not satisfied (figure |4p corre- 
sponds to this case) , any horizontal line e = const crosses 
some green region, i.e., there are always waves that can 
penetrate through the barrier. 



6 



What is important for potential applications is that the 
parameters of the barrier can be easily controlled by the 
applied voltage. Depending on the voltage, this barrier 
can be either semitransparent (Fig. \%jp) or opaque [with 
bound states made of surface waves (Fig. SH)]- 



IV. DOUBLE BARRIERS CAN PRODUCE 
WAVEGUIDES 

Under certain conditions, two barriers separated by a 
distance d form a waveguide. Two barriers divide the 
infinite graphene sheet into tree domains: left and right 
infinite half-planes, and the strip (of width d) bounded 
by the barriers. The potentials (electric and magnetic) 
are constant within each domain. The values of the po- 
tentials are denoted by the sub-indices £, r, and c, respec- 
tively, for left, right, and middle. Each domain is repre- 
sented by its own division of the (k±, e)-plane in sectors 
with propagating and non-propagating (in the given do- 
main) waves, as it was described in the previous section. 
The situation is similar to those presented in Fig. Hti,b, 
with only one difference: now there are tree Dirac points 
(—Aa, u a ) (a = £,c,r) with "cones" that divide the 
plane on several sectors. Hereafter, the yellow-colored 
sectors correspond to the central (c) graphene domain, 
the blue-colored and violet-colored sectors correspond to 
the right (r) and left (I) graphene domains, respectively. 
All overlaps of the sectors will be marked by dark blue 
(£— r overlap), light green for the i — c overlap, green for 
t—c overlap, and dark green for r-c-l overlap. 

The waveguide eigenmodes are non-propagating, 
evanescent waves in the left and right infinite graphene 
half-planes. This means that the corresponding points in 
the (k±, e)-plane are located outside of the colored sec- 
tors which are generated by the r and I Dirac points, i.e., 
in either yellow or white regions of the (k_l, e) -plane. 

To simplify the presentation, we will mainly consider 
symmetric cases, \ug\ — \u r \, and \Ai\ — \A r \ and assume, 
without loss of generality, that u c = A c = 0. 



A. Equal scalar and vector potentials: ui = u r and 

At = Ar 

We first consider the case when ug = u r and Ag — 
A r . Setting ui — u r — U2 and At = A r = A2 (the 
blue-colored sectors correspond to the r and t equivalent 
graphene domains), and setting u c = u% = and A c — 
Ai = for the central domain, one can now use Fig. 2k,b 
for describing the double-barrier structure. 

When the inequality (j2"Tj) is not valid (as in Fig. 0b), 
there are waveguide modes which are similar to the 
modes in the usual dielectric waveguides. The spectra 
of these modes are shown by the red lines in Fig. (S£i,b 
for different values of the distance d between the barriers 
(waveguide walls). 



As in a dielectric waveguide, the total internal reflec- 
tion (TIR) phenomenon is responsible on the wave con- 
finement in the central graphene domain. However, in 
contrast to usual dielectric waveguides, in a certain en- 
ergy range there are two separate regions of k± where 
TIR occurs. Decreasing the distance d between the bar- 
riers shifts the spectrum to higher-energies, with the 
exception of the lowest mode, which crosses the point 
{—At, ui) for whatever small values of d (Fig. ^p). 



(a) e 




FIG. 5: (color online) Waveguide spectrum e(ftx) (red lines) 
when the inequality (1211) is not valid, (a) for large spacing 
d between the barriers; (b) for small spacing d between the 
barriers. 

When the inequality (|21l) is valid, there are confined 
waveguide modes in the yellow regions in Fig. Ht 1 - In addi- 
tion to these "ordinary" modes, there are two "extraor- 
dinary" modes in the white region between two Dirac 
points where the surface wave exists. In contrast to the 
ordinary modes that appear due to the TIR, the extraor- 
dinary modes are formed by two coupled surface waves 
propagating along the barriers. 

In Fig. [Bk, the blue 2D "cone" on the right represents 
the overlap of two 2D "cones" corresponding to the left 
and right graphene semiplanes. The isolated red solid 
line originated at the right Dirac point actually repre- 
sents two dispersion curves, which in this instance are 
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indiscernible. When these two nearly-overlapping curves 
approach the left Dirac point, they separate moving in 
opposite directions: one towards positive e, the other 
towards negative e. In other words, the two coupled 
surface waves propagating along the barriers (red solid 
curve emanating from the right Dirac point) transform 
smoothly into the standard ordinary modes, when ap- 
proaching the left Dirac point. Examples of waveguide 
spectra are shown in Fig. [6£i,b for different values of the 
distance d. 

The particular case up, = u r = has been considered 
in [ill IH, l32l . l33t , a nd the case Ap = A r = has been 
considered in [34], [35[ . 



(a) £ 




(b) £ 




FIG. 6: (color online) Waveguide spectra (red lines) when the 
inequality (|2ip is valid, (a) large spacing d between the bar- 
riers. The isolated red line in (a) involves two extraordinary 
modes, corresponding to edge states, which merge in the iso- 
lated red line because these two curves are very close to each 
other. The difference between these nearly-overlapping modes 
is visible only when the modes are transformed into ordinary 
ones, (b) small spacing d between the barriers. Note that the 
spacing between the energy levels increases when the spacing 
decreases. 

The eigenspectrum of both ordinary and extraordinary 
modes is defined by the requirement that a round trip 



phase 

cfi = 2Re(K||)e? + (f>e + (j> r , 

(<fip and <fi r are the phase shifts at the reflection from 
the left and right barriers) is either equal to zero (ex- 
traordinary evanescent modes) or to a multiple of 27m, 
n = 1,2,... (ordinary propagating modes). 

There is an important difference between the ordinary 
and extraordinary modes: when the distance between 
the barriers is large enough, the extraordinary modes 
are unidirectional waves (see Fig. [5^). This means that 
all the waves (in the extraordinary mode existence re- 
gion) propagate in one direction, making this mode re- 
sistant to backs cattering, and therefor robust against y- 
dependent disorder. Indeed, any rather smooth obsta- 
cle in the waveguide cannot produce counter-propagating 
waves. As to the fluctuations of the height of the poten- 
tials (in a;-direction) , in principle, they can affect dra- 
matically the propagation across the barriers, but have 
practically no influence on the waveguide eigenmodes. 
Note that the properties of the waveguide eigenmodes 
are easily controlled by the electrostatic potentials. 

The closeness of two extraordinary waves spectra (cen- 
tral isolated red line in Fig. [5^) stems from the chosen 
symmetry of the waveguide potentials up — u r . When 
this symmetry is broken, up ^ u r , the difference be- 
tween the spectra is clearly visible (two central red lines 
in Fig [7^). When the potentials up and u r have opposite 
signs, surface waves propagate along the barriers in oppo- 
site directions. As a result, two extraordinary modes not 
only have separated spectra, but propagate in opposite 
directions, as it is shown in Fig [7)d. The backscatter- 
ing is also suppressed in this case because the counter- 
propagating waves are localized mainly near the corre- 
sponding barriers. 

B. Equal scalar potentials and antiparallel vector 
potentials: up — u r and At = — A r 

In this subsection, we will concentrate on the case when 
Eq. (|2"Tj) is valid. The instance when the inverse of the 
(|2ip inequality is valid is similar to one considered in 
subsection IV A, just with a more complex division of 
the {k±, e)-plane in regions. 

The condition up/Ap — —u r /A r means that the sur- 
face waves along the left and right barriers propagate 
in opposite directions making the waveguide eigenmodes 
spectra symmetric, as shown in Fig. [5£i,b. 

In contrast to the case of equal potentials (subsection 
IV. b) the spectra of ordinary and extraordinary modes 
are separated: there is no transformation of one kind 
of mode into another, as it is seen in Fig. [S] and Fig. [7J 
Note that the coupling between two surface waves, which 
produces the extraordinary mode does not split the mode 
energy, as it happens when up = u r and At = A r . 

There is an energy range where only extraordinary 
modes form the charge flux along the waveguide. De- 



8 





FIG. 7: (color online) Waveguide spectra for the cases dis- 
cussed in the subsection IV A. (a) for ue 7^ u r , sgn(wf) = 
sgn(« r ); (b) for ui = —u r . The difference between the 
two extraordinary modes is distinctly visible. For the case 
sgn(tif) = — sgn(tt r ), the modes in (b) have opposite-directed 
group velocities. 



FIG. 8: (color online) Waveguide spectra when the scalar ut = 
u r and vector At = —A r potentials satisfy these relations, (a) 
large distance d between the barriers; and (b) small distance 
between the barriers. When the distance d decreases, the 
ordinary modes are pushed out from the yellow regions and 
only the extraordinary mode remains in a broad energy gap. 



spite the fact that the extraordinary modes are not uni- 
directional, the backscattering is also suppressed in this 
case. The reason for this is the spatial separation of re- 
gions with opposite direction of the flux: the wave func- 
tions that correspond to opposite directions of flux are 
localized near different (left or right) barriers. This is 
illustrated in Fig. [9j where the normalized flux densities 
for extraordinary eigenmodes with opposite wave vectors 
k±i = —k±2 (opposite total fluxes) are shown. 

Decrease of the spacing, d — > 0, broadens the energy 
gap where only extraordinary mode exists (see Fig. [9Jd) . 
The dispersion curve is then flattened and tends to the 
line s(kx.) = up = u r . 



C. Scalar potentials of opposite signs and 
antiparallel vector potentials: U( = —u r and Ai = —A r 

When up — —u r , At = —A r and the inequality (p?Tj) 
holds, both barriers support surface waves with the same 



direction of the charge flux. Accordingly, the flux which 
is associated with the extraordinary waves is unidirec- 
tional, i.e., is independent on sign of k±. Examples of 
waveguide spectra are shown in Fig. [TUa.b. When the 
distance d is small enough, only extraordinary waves are 
confined by the waveguide. 

Comparing Figs. [TUb andllOb. one can see that, unlike 
the ordinary modes, the spectrum of the extraordinary 
modes is independent of the distance d between barri- 
ers. There is a cut-off energy £0 for ordinary modes, 
i.e., the minimal modulus |e| of the energy, for which the 
condition </> = 2ir is satisfied (the cut-off energy is the 
analog of the cut-off frequency of conventional waveg- 
uides). The cut-off energy increases when the distance 
d decreases, pushing out the ordinary modes from their 
existence region (yellow regions in Fig. [TUk.b). It is in- 
teresting to note that in the graphene waveguide with 
variable width d(jj), the extraordinary mode can pen- 
etrate through an arbitrary narrow part of the waveg- 
uide, whereas all ordinary modes are reflected from it. 
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FIG. 9: (color online) Normalized current density of the 
extraordinary modes (counter-propagating surface waves). 
Solid line: k± > 0; dashed line: k± < 0. The vertical lines in- 
dicate the positions of the barriers. The counter-propagating 
currents are indeed spatially separated. 

A similar effect (penetration of the electromagnetic wave 
through the waveguide waist) is typical to waveguides 
filled with a metamatcrial with a near-zero dielectric per- 
mittivity [36j . A rather narrow waveguide is a single- 
mode waveguide, although its transport properties differ 
strongly from the usual single-mode waveguides. The 
transmission through this waveguide (in the frame of the 
model we use) is always perfect, irrespectively of smooth 
irregularities present on the waveguide. As in the sin- 
gle barrier case, the reason is the absence of counter- 
propagating waves that makes reflection impossible. 

V. CONCLUSIONS 

We have shown that crossed magnetic and electric 
fields applied to a narrow strip (electromagnetic barrier) 
on a graphene sheet can form a unidirectional conduct- 
ing channel (quantum wire) whose properties are easily 
tunable by the voltage applied across the strip. The 
eigenmode of this channel is characterized by a linear 
dispersion and represents a one-way propagating wave. 
This unique property prohibits backs cattering and there- 
fore makes the mode resistant to the scattering by impu- 
rities. The classical analogy of this mode is the drift 
of a charged particle in crossed electric and magnetic 
fields. The transverse-localized mode exists and prop- 
agates longitudinally along the barrier with the drift ve- 
locity v y — Vd = cE x /H z if and only if this velocity Vd is 
smaller than the Fermi velocity Up- i.e., Vd < vp. While 
one barrier forms a wire, two such barriers produce a 
waveguide whose eigenfunctions consist of set of ordinary 
and extraordinary waves. The ordinary waves are char- 



acterized by the quantization of their transverse wave 
numbers, while the extraordinary waves are formed by 

(a) e 




(b) s 




FIG. 10: (color online) Waveguide spectra for the case de- 
scribed in subsection IV.C, with opposite scalar and vector 
potentials, (a) large spacing between the barriers; (b) small 
spacing between the barriers. 

two coupled surface waves propagating along the waveg- 
uide walls (barriers). The rather narrow waveguide has 
only extraordinary eigenmode. Depending on the param- 
eters of the walls this extraordinary eigenmode can be ei- 
ther uni- or bi-directional. In the bi-directional case the 
regions with opposite directions of current flows are spa- 
tially separated, preventing backscattering and making 
even the bi-directional mode resistant against the scat- 
tering by impurities. 
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